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The discrepancy between the amplitudes of matter fluctuations inferred from Sunyaev-Zel’dovich 
(SZ) cluster number counts, the primary temperature, and the polarization anisotropies of the 
cosmic microwave background (CMB) measured by the Planck satellite can be reconciled if the 
local universe is embedded in an under-dense region as shown by Lee, 2014. Here using a simple 
void model assuming the open Friedmann-Robertson-Walker geometry and a Markov Chain Monte 
Carlo technique, we investigate how deep the local under-dense region needs to be to resolve this 
discrepancy. Such local void, if exists, predicts the local Hubble parameter value that is different 
from the global Hubble constant. We derive the posterior distribution of the local Hubble parameter 
from a joint fitting of the Planck CMB data and SZ cluster number counts assuming the simple void 
model. We show that the predicted local Hubble parameter value of iLioc = 70.1±0.34 km s“^Mpc“^ 
is in better agreement with direct local Hubble parameter measurements, indicating that the local 
void model may provide a consistent solution to the cluster number counts and Hubble parameter 
discrepancies. 

PACS numbers: 98.70.Vc, 95.30.-k, 98.80.Es 


I. INTRODUCTION 


The standard model of cosmology, also referred to the 
ACDM model, is now well established based on a num¬ 
ber of cosmological observations. The simplest ACDM 
model contains six free cosmological parameters, each of 
which has been determined with an accuracy of a few 
percent through the measurement of Cosmic Microwave 
Background (CMB) anisotropies by the Planck satellite 
[l| . The CMB anisotropies mainly manifest the nature of 
density fluctuations in the universe at the recombination 
era at 2 1100. It is therefore of great importance to 

confront the standard model with observations of the lo¬ 
cal universe {z ss 0) to confirm that the standard model 
is correct throughout the history of the universe (see sec¬ 
tion 5.6 of 0 and references therein). 

Some authors, however, have reported cosmological 
measurements in the local universe that conflict with the 
CMB result. For example, it has been claimed that the 
amplitude of density fluctuations measured in the local 
universe is systematically low compared with the predic¬ 
tion from Planck’s CMB data (e.g., [1-Q)- In particular, 
the number of massive clusters derived from the Sunyaev- 
Zel’dovich (SZ) effect by Planck is about half that ex¬ 
pected from the CMB anisotropies M- This SZ result 
from Planck is consistent with the SZ-selected cluster 
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number counts from the Atacama Cosmology Telescope 
m and the South Pole Telescope surveys [T^. Perhaps 
related to this discrepancy, the local Hubble parameter 
measured through the traditional distance ladder tends 
to show a larger value (e.g., Hq = 73.8±2.4kms“^Mpc“^ 
by m) than that inferred from the Planck result {Hq ^ 
67.3 ± 1.0 kms“^Mpc“^). There have been many pro¬ 
posals to resolve these discrepancies, including beyond- 
the-standard models, e^., massive neutrinos [j, and 
decaying dark matter |14| , and re-considerations of astro¬ 
physics and calibration issues, e.g., the halo mass func¬ 
tion [IM3 , and mass calibration [isl-l^. 


Here, we reinvestigate the idea of Lee 2014 [23|, where 
it was shown that the discrepancy can be resolved if we 
reside in a local under-dense region. The idea is partic¬ 
ularly interesting because the local under-dense region 
may explain the observed discrepancy in the Hubble pa¬ 
rameter as well [ 2 ^. In our analysis, we employ a sim¬ 
ple Friedmann-Robertson-Walker (FRW) model with an 
open geometry for the under-dense region and investi¬ 
gated how deep the under-dense region is in order to 
resolve the discrepancy using the Markov Chain Monte 
Carlo (MCMC) technique. We show that the joint fitting 
of the CMB and SZ data set predicts the local Hubble 
parameter value. The posterior distribution of the lo¬ 
cal Hubble parameter from the MCMC analysis is then 
confronted with direct local Hubble parameter measure¬ 
ments to check the consistency of the local void picture. 


The paper is organized as follows. In the next sec¬ 
tion we describe our simple void model and our method 
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of analysis. We present our results and discussions in 
Section III. The final section is devoted to our conclu¬ 
sion. The normalized Hubble parameter h is defined by 
h = Ho/{100 kms“^Mpc“^). 


II. METHOD 
A. Local void model 


Suppose that we live in a local (under-)dense region in 
a background flat Friedmann-Robertson-Walker (FRW) 
universe. If the local region has a constant density pro¬ 
file, i.e., a top-hat profile, the evolution of the region 
is approximated by a slightly different FRW universe. 
The local universe and the background universe may 
be characterized by cosmological parameters:(iLioc, flioc) 
and (iJbgjf^bg), respectively. The local cosmological pa¬ 
rameters (FIioc, f^ioc) can differ in principle from those in 
the global FRW universe. Here, for simplicity, let us as¬ 
sume that abundance ratios of the energy components in 
the universe are same in the global and the local under- 
dense regions by demanding that 


f^a.loc^loc — 


( 1 ) 


where the subscript a = (v, c, b, 7 , n) stands for the vac¬ 
uum energy (v), CDM density (c), baryon density (b), 
photon density ( 7 ) and neutrino density (n). Note that 
we have assumed that the curvature parameter in the 
background FRW model is equal to zero, i.e., Hx.bg = 0. 
The local curvature parameter, Hk, 1 oc is determined by 


f^K.loc — 1 / . ) (2) 

i 

so that the difference between H\oc and iJbg is compen¬ 
sated. Then, we find Hk, 1 oc > 0 (< 0) for H\oc > H^g 
{Hioc < Hhg)- It should be noted that the present time 
must be fixed so that the present CMB temperature is 
2.725K. Therefore, our assumption o implies that, for 
the present uniform density time slice, the CMB temper¬ 
ature we observe and that in the background universe are 
equally given by 2.725K. Equation © also implies that 
the local void considered here is in the growing mode (see 
e.g., Eq.(34) in [H]). 

Let us introduce a parameter /loc as follows: 


Hioc = fiocHhg ■ (3) 

We are particularly interested in the case /loc > 1 to ex¬ 
plain the discrepancy in the measured values of the Hub¬ 
ble parameter. In this case, in order to change the present 
time slice from the uniform density slice to the uniform 
Hubble slice given hy H = Hioc > .ffbg, we have to go 
back in time and take the past slice in the background 
universe. Consequently, in the background universe, the 
density on this time slice is larger than that on the uni¬ 
form density time slice. This means that the density 


profile has a void structure in the uniform Hubble slice. 
This observation is consistent with the result reported in 
Ref. [ 2 ^, which states that the local under-dense region 

may lead to a larger value of the local Hubble parameter. 
1 

In this paper, we calculate the cluster abundance based 
on the local cosmological parameters: (iJioc, f^ioc), which 
are fixed by (i7bg,f2bg) and /loc through Eqs. o and 
([3|). In other words, we assume that the local void cov¬ 
ers up the entire region where the Planck clusters were 
found, namely, z < 1 , or in terms of the radial comoving 
distance, < 3000 Mpc. 

The effect of the parameter fioc on the cluster abun¬ 
dance is shown in Fig. [T] The decrease in the clus¬ 
ter abundance with increasing /io<. mainly comes from a 
smaller amplitude of density fluctuations at correspond¬ 
ing scales due to a slower growth rate in the open FRW 
model than that in the flat FRW one. Following the con¬ 
vention, we denote the fluctuation amplitude by ctsjoc, 
the root-mean-squared of linear fluctuations within a top- 
hat sphere of 8h~^ Mpc radius. Assuming that primor¬ 
dial density fluctuations have the same initial power spec¬ 
trum both in the under-dense region and the surrounding 
background universe, the parameter crg^oc is derived from 
the other cosmological parameters. 



redshift 


^ The correlation between the larger value of the local Hubble 
parameter and the local under-dense region can be also under¬ 
stood by considering conventional perturbation theory. Usually, 
the initial condition for the Hubble flow of the top-hat region is 
set assuming that the density perturbation grows following lin¬ 
earized perturbation theory; in the radiation dominated era, the 
density contrast of region <5 grows as <5 oc in the synchronous 
gauge, where a is the cosmic scale factor. In this case, the re¬ 
lationship between the Hubble parameters of the top-hat region 
Hioc a-ud the background universe i7bg is given by 1261 

2 

HloCji Hlog^i , 

where the subscript i denotes a specific initial time. If the top- 
hat region is under-dense, i.e., —1 < <5 < 0, then H\oc,i > H\^^ i 
and the inequality H\oc{t) > Ilbg(f) always holds true. 
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FIG. 1: SZ cluster number counts predicted in the local void 
model are compared with the observed number counts from 
the Planck 2013 result 0- The parameter / = /loc defined in 
Eq. m describes the difference between the local and global 
Hubble parameters. The parameter h here denotes the nor¬ 
malized local Hubble parameter value corresponding to each 

/loc 


For simplicity, we assume that the angular diameter 
distance to the last scattering surface remains unchanged. 
This can be justified if the density profile in the high red- 
shift universe, say in 2 < z < 1000, is modified to com¬ 
pensate for the stretch in distance by the local under- 
dense region. In this case, the Hubble parameter Ffbg is 
not directly observed but inferred from the CMB mea¬ 
surements. On the other hand, i?ioc can directly be ob¬ 
served using local distance measurements. 


z 

Mnim(z) [Mq\ 

0.05 

3.12el4 

0.15 

6.15el4 

0.25 

8.20el4 

0.35 

9.66el4 

0.45 

1.06el5 

0.55 

1.13el5 

0.65 

1.17el5 

0.75 

1.20el5 

0.85 

1.21el5 

0.95 

1.21el5 


TABLE I: Mean minimum masses of the Planck 2013 SZ clus¬ 
ter sample as a function of redshift, which are taken from 0. 
These masses are defined at 50% completeness on average for 
the unmasked sky. 


B. Mass function and cluster abundance 


To predict the abundance of SZ clusters, we adopt 
the fittin g fo rm of the mass function drihaio/dfl/ pre¬ 
sented in [^. We adopt the parameters Hcrit = 500 and 
£*ciuster = Dcrit/^Miz), where Hm( 2 :) is the total matter 
density over the critical density at redshift x. The cluster 
abundance of mass M at redshift z is given by 


dN, , , r ^dV{z) 

= "xlM) —(M,.)—, (4) 


where the co-moving volume V(z) is given by 


dV(z) 

dz 


= 47r(l -I- z)^ 


dAizf 

H{z) 


(5) 


and the sky fraction that is covered by the Planck SZ 
cluster counts is fixed to /sky = 0.65. The function x(M) 
represents the survey completeness as a function of halo 


mass M that takes account of the fact that some frac¬ 
tion of clusters falls out from detection due to the survey 
strategy, detection limit, and so on. We simply write this 
function as [2^ 


X{M) = / dM'P{M'\M) , (6) 

where the function P{M'\M) describes the probability 
that the mass is estimated as M' for a cluster with a true 
mass M, which for example originates from the scatter in 
the scaling relation between halo masses and SZ signals. 
In this analysis, we assumed that it follows a log-normal 
distribution with variance ^ as 


P{M'\M) 


M'-i 


exp 


^O'ln M 


(7) 


We fix tTinM = 0.2, which roughly reproduces the theo¬ 
retical curve for the number counts of Planck SZ clusters 
presented in [0|, as shown by the black curve in Fig. [T] 
The minimum masses obtained from the Planck SZ sur¬ 
vey are listed in Table HI 

A caveat is that, while we have used the Tinker mass 
function to compare our results with those of the Planck 
fiducial analysis, the fitting function is not very well 
tested in an open ACDM model which is of our main in¬ 
terest in this paper. We expect that the fitting function 
works fine even in non-flat universes as long as important 
ingredients of the fitting function such as the matter den¬ 
sity, linear power spectrum, and linear growth rate are 
properly included, but this assumption should be care¬ 
fully validated by e.g., Wbody simulation. 


C. MCMC 

In our MCMC analysis we use both the Planck CMB 
data and the SZ cluster number counts. Specifically, we 
used a flat ACDM model when fitting to the CMB data 
with the six standard parameters, while we used an open 
ACDM model with an extra parameter /loc to ht to the 
SZ cluster data; the other parameters were derived from 
the relations in Eqs. ®-©- We use a modified version 
of CosmoMC [ 2 ^ to explore the likelihoods. For the like¬ 
lihood of the cluster number counts, we assume that the 
number of the cluster obeys Poisson statistics [s^ , which 
is a good approximation for massive clusters as consid¬ 
ered here 13 ill. 


III. RESULTS AND DISCUSSION 

First, we show the posterior distribution for /loc from 
the joint fitting of the CMB and SZ number counts in 
Fig. [5] As shown in the Figure, the parameter fy^c tends 
to be /loc > 1 indicating that the model with a local 
under-dense region fits better with the data set. The 















4 


result /loc = 1-03 ± 0.01245 indicates that the standard 
ACDM model is excluded at 2.4tT level. We find the im¬ 
provement of to be Ax^ 8.44 with one additional 
parameter, which is, in terms of Akaike’s information cri¬ 
terion, AIC ~ 6.44. This indicates that the void model 
is preferable to the ACDM model for the combined CMB 
and SZ data set. The two dimensional constraint on the 
/loc — cs.ioc plane is shown in the left panel of Fig. |31 It is 
evident that the trend for /loc > 1 is required in order to 
realize the smaller ag as inferred from the measurements 
of the local universe. 

An interesting relation is found in the right panel of 
Fig 131 which shows a positive correlation between crg^ioc 
and i?ioc- This correlation is non-trivial because a larger 
floe leads a smaller crgjoc but a larger Hioc if the other 
cosmological parameters are fixed. We find that this cor¬ 
relation is attributed to the well-known crg-Dm degener¬ 
acy in the cluster abundance. In order to produce the 
same cluster abundance, a larger crg^ioc requires a smaller 
f^m,ioc 5 hence a larger Hioc because 
is tightly constrained from the CMB. 




*^8,loc ^8,loc 

FIG. 3: Posterior probability distribution on the /loc — cg^ioc 
{left) and Hioc — crg^ioc (right) planes. 


parameter is predicted by the joint fitting of the CMB 
and SZ cluster number counts as Hioc — lOOhioc = 
70.1 ± 0.34 kms“^Mpc“^. This value is in better agree¬ 
ment with the local Hubble measurements by Riess et al. 
01 and interestingly, is in excellent agreement with the 
value Hq = 70.6 ± 3.3 kms“^Mpc“^ recently derived by 
Efstathiou usiM the updated geometric maser distance 
to NGC4258 [3|. 



62 64 66 68 70 72 74 64 66 68 70 72 74 76 78 


Hq [km/s/Mpc] 


Hq [km/s/Mpe] 


FIG. 2: One-dimensional posterior probability distribution 
of the parameter /loc (see Eq. [3| from the joint fitting of 
the CMB and SZ cluster number counts. The dotted line 
(/loc = 1) denotes the standard cosmological model without 
any local over- or under-dense region. 

A key feature of this local void model is that the lo¬ 
cal Hubble parameter defined in Eq. is automatically 
adjusted as a result of fitting to the observed SZ clus¬ 
ter number counts. Figure U shows the posterior prob¬ 
ability distributions of the Hubble parameters with and 
without the void parameter /loc, together with the local 
Hubble parameter measurements from Riess et al. [l^ 
and Efstathiou (s^ . The main results are summarized 
as follows. By keeping the void parameter at the stan¬ 
dard ACDM value (/loc = I), the Planck CMB gives 
Hq = 67.3 ± 1.2 kms“^Mpc“^ whilst the CMB and SZ 
combination gives Hq = 70.7 ±0.3 kms“^Mpc“^, indi¬ 
cating the disparity between these two data sets. We 
find that in the current void model, the local Hubble 


FIG. 4: (Left) Probability distribution functions for the Hub¬ 
ble parameter from Planck (black solid line), and Hbg (blue 
solid) and Hioc (blue dashed) from the combined CMB and 
SZ data set. (Right) Constraints on the Hubble parameter are 
compared with the direct local estimates for the Hubble con¬ 
stant by Riess et al. (Rll), and by Efstathiou using NGC 4258 
(E14) and plus other distance anchors (E14-I-MW-I-LMC). 
The black point with an error bar shows the constraint from 
Planck CMB. In the local void model, Hubble parameters are 
separated into local (Hioc) and global (Hhf) parameters, as 
shown in the figure. 

There are several other possibilities to resolve the dis¬ 
crepancy. For example it may be possible that the mass 
bias factor is significantly lower than the inferred value, 
i.e., masses of clusters observed by Planck SZ is in fact 
larger than those inferred by the scaling relation, and 
therefore, they are less abundant in the universe [l^. 
Since the publishing of their results in 2013, the Planck 
collaboration has updated their bias calibration using 
cosmic shear and CMB lensing measurements, and their 
2015 estimates give similar mass bias values to that in 
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the baseline model assumed in 2013 [s^. The results of 
cluster number counts in 2013 and 2015 are consistent 
with each other if the same mass bias is assumed. In¬ 
dependent analysis also shows that the mass bias cannot 
fully resolve the discrepancy Recently, an im¬ 

portant update has been made in Battaglia et al. (^ . 
They found that an Eddington bias correction, if it is 
applied to the weak lensing mass calibration analysis of 
the Planck SZ clusters, brings the mass bias factors mea¬ 
sured in von der Linden et al. [l^ and Hoekstra et al. 

closer to the values that are small enough to explain 
both the Planck CMB and SZ results. Further analysis 
of this systematic uncertainty is of primary importance 
for understanding the origin of the discrepancy. 

From a cosmological point of view, while massive neu¬ 
trinos may be an interesting possibility an, a larger 
neutrino mass generally demands a smaller Hubble pa¬ 
rameter, which makes the discrepancy between the global 
and local Hubble parameters even worse [s^l . Other pos¬ 
sibilities include modification of the mass function due to 
baryonic effects [TB - [l7l| . and decaying dark matter [l3| . 
Non-Gaussianity may also explain the discrepancy, al¬ 
though the Planck bi-spectrum already put severe con¬ 
straints [ 3 ^, making this explanation difficult as long as 
the non-Gaussianity is scale-independent [^ . 

Our current study needs to be advanced in several di¬ 
rections. The large bulk flow associated with the local 
void generates secondary GMB anisotropies through the 
kinetic Sunyaev-Zeldovich effect (kSZ). Here we give a 
rough estimate of the kSZ effect arising from the lo¬ 
cal void following Moss et al. [i^. Suppose that the 
dipole induced by the radial bulk velocity is the domi¬ 
nant anisotropy seen by the scatterer (i.e., clusters). We 
estimate the amplitude of this dipole from the difference 
in redshift between incoming and outgoing CMB photons 
felt for the same period in conformal time. The dipole 
anisotropy at clusters is then given by 


^ ~ 1 - 
T a{ri^,) 


( 8 ) 


where 0 ( 77 *) and R)??*) denote the scale factors of the 
background and void regions when the incoming photon 
enters into the void, and follow the Friedmann equations 
in flat and open ACDM models, respectively. For our 
void model that gives H\oc = 70, we obtain AT/T « 
7 X 10-4 for a(? 7 *) w 0.01 and AT/T r; 9 x lO’^ 
for a(? 7 ») w 0.1. Therefore, a fairly large void seems 
necessary to avoid the constraint from peculiar veloci¬ 


ties measured by the Planck collaboration, which gives 
AT/T < 6.4 X 10-4 (2cr) [dlj. These numbers above 
are order of magnitude estimates; further detailed in¬ 
vestigation would be needed for a more quantitative con¬ 
straint by using concrete models of spherically symmetric 
inhomogeneous universes, e.g^ generalized A-Lemaitre- 
Tolman-Bondi dust models [4^ - 1^ . The effect on other 
probes such as baryon acoustic oscillation and cosmic 
shear measurements should also be investigated in or¬ 
der to see whether the local void is a viable cosmological 
model. 

This study is based on the Planck 2013 results. In¬ 
terestingly, the discrepancy persists in the Planck 2015 
results and therefore our results should still be as 
valid for the latest data set. 


IV. CONCLUSION 


In this paper, we have considered a simple FRW model 
with an open geometry to investigate the magnitude of 
the effect of the local under-dense region on the number 
of SZ clusters. We have shown that a local under-dense 
region may be responsible for the discrepancy between 
cosmological parameters from SZ cluster number counts 
and the CMB. If we attribute the discrepancy to the 
local void, the local Hubble parameter is predicted as 
Hq = 70.1±0.34 kms-4Mpc-4, which is in better agree¬ 
ment with direct local Hubble parameter measurements. 
Our result indicates that the local void model may serve 
as a new concordance model that explains the CMB, SZ 
cluster number counts, and local Hubble parameter mea¬ 
surements simultaneously. 
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